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^ Abstract 

We seek to gain insight into the nature of the determinantal moments | \p\ ) of 

O generic (9-dimensional) two-rebit and (15-dimensional) two-qubit systems (p), PT denoting par- 

tial transpose. Such information-as it has proved to be in the Hilbert-Schmidt counterpart-should 

be useful, employing probability-distribution reconstruction (inverse) procedures, in obtaining En- 
s' 

Q_i res 2x2 separability probabilities. The (regularizing) strategy we first adopt is to plot the ratio 

^ of numerically-generated (Ginibre ensemble) estimates of the Bures moments to the corresponding 

^ (apparently) exactly-known Hilbert-Schmidt moments (J. Phys. A, 45, 095305 [2012]). Then, 

through a combination of symbolic and numerical computations, we obtain strong evidence as to 

the exact values (and underlying patterns) of certain Bures determinantal moments. In partic- 
le 

Q\ ular, the first moment (average) of (where \p^'^\ G ^fel) two-qubit systems 
(N 

T-H is, remarkably, — jgg = ~2 ^ ~ —0.00390625. The analogous value for the two-rebit systems is 

O -silo = ~ 2i^Sx7 - -0.00309594. While in the Hilbert-Schmidt case is the ratio 

I of 3n-degree polynomials in k, it appears to be the ratio of 5n-degree polynomials in k in the Bures 



, case. 

X 
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I. INTRODUCTION 



Zyczkowski, Horodecki, Sanpera and Lewenstein raised and investigated-in two-qubit and 
qubit-qutrit settings-the clearly non-trivial, fundamental quantum-information-theoretic 
question of what proportion of generic bipartite quantum states is separable [1]. Over 
the past several years, we had-employing Bures pi Ej and Hilbert-Schmidt |3l H] measures 
on the space of such states-examined this question using several distinct approaches-but 
with no definitive resolutions [5HT2] . The problem appears to be particularly challenging, in 
large part, due to the high-dimensionality of the relevant spaces-15- and 35-dimensions in 
the generic two-qubit and qubit-qutrit cases, respectively. 

Recently, we have reported substantial advances — in the two-qubit (2 x 2) case, employ- 
ing Hilbert-Schmidt (HS) measure-using a (somewhat indirect, but nevertheless effective) 
multi-stage approach based on, first, obtaining formulas for the HS moments of 4 x 4 density 
matrices (p) and of their partial transposes [T3HT5]. (For example, the evidence derived- 
though not yet fully rigorous-is compelling that the Hilbert-Schmidt separability probabili- 
ties of generic [9-dimensional] two-rebit states is ||, of [15-dimensional] two-qubit states is 
|j, and of [27-dimensional] two-quater(nionic)bit states is ^.) It is the objective of this 
study to investigate and evaluate the possibility of fully adopting this apparently successful 
approach to the Bures well. We are able to report some interesting progress in this 

direction, but still fall considerably short of such a long-range objective. We will now briefly 
discuss these recent (Hilbert-Schmidt-based) advances, that motivate our present quest to, 
additionally, find Bures analogues. 

In [131 App. D.6], formulas (not yet rigorously established) were developed by Dunkl for 
the HS bivariate product determinantal moments (jp^^\" \p\^^ , fc, n = 0, 1, 2, 3, . . . ,-PT 
denoting the partial transpose. Here, a is a Dyson-index-like parameter, with, in particular, 
a = 1 denoting the two-qubit systems and a = ^, the two-rebit systems. These formulas 
then served as input (using the induced univariate moments (|p^"^|")^^ ^-^id (Ip"^"^!" \p\"^)^^ 
[T3| E]) to a (Legendre-polynomial-based) probability-distribution reconstruction proce- 
dure [16]. This provided estimates (and eventually-ma a hypergeometric-based expression 
generated-exact results) for Hilbert-Schmidt 2x2 separability probabilities as a function of 
a. 

In sec. [IT] we will compare-graphically-Bures and Hilbert-Schmidt 2x2 determinantal 
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FIG. 1: Ratio of numerically-generated Bures moments (|p-'^^|")^^^^*^^_,-based upon 59,800,000 
random realizations-to corresponding exact Hilbert-Schmidt moments 



moments for generic rebit and qubit systems, while in sec. Ill we will begin-what appears to 



be a laborious process-of obtaining exact results pertaining to Bures determinantal moments. 



II. BURES/HILBERT-SCHMIDT DETERMINANTAL MOMENT RATIO PLOTS 

In the present study, we first employed the Ginibre ensemble procedures detailed in [TTj 
to generate 59,800,000 (9-dimensional) generic two-rebit and (15-dimensional) generic two- 
qubit density matrices distributed according to Bures measure. Then, using high-precision 
arithmetic, we obtained estimates of the first 500 moments of the determinant |p'^"^|-the 
nonnegativity of this determinant being necessary and sufficient for separability in the 2x2 
case [IB]. 

In Figs.[T]and[2]we plot the ratios of these estimates to the exactly-known Hilbert-Schmidt 
moments using one of the formulas presented in [131 App. D.6], that is 

,| pT^n.HS ^ n\ (a + l)„(2a + 1)„ 
' 2-(3a + i)„(6a + i),^ 

(_2n-l-5«)„« ^ / _n_2,_n_i,_^,« + i,2a + l \ 

24" (3a+ 1)^ (6a+ ^ ^ \1 - n,n + 2 + 5a,l - n - a, ^ - n - a J' 

(Because of the occurrence of the term 1 — in the denominator, it is necessary to replace 
the 5F4-sum by 1 to obtain the correct value when n = 1.) This graphic approach might be 
regarded as a form of regularization of our numerical procedures, with the resulting ratios 
plotted in Figs. [T] and [2] being of far less extreme magnitudes than the numerical estimates 
of the Bures moments themselves. 
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FIG. 2: Ratio of numerically-generated Bures moments (|p^^|")^_"^u^jts~^^®^'^ upon 59,800,000 
random realizations-to corresponding exact Hilbert-Schmidt moments 



Since the lower-order Bures moments dominate the Hilbert-Schmidt ones in absolute 
value, it would appear that it would-attractively-take fewer Bures moments (if known ex- 
actly) than those several thousand Hilbert-Schmidt ones employed in [131 Ej to obtain 
comparable accuracies in the probability-distribution reconstruction procedure [16]. 

Further, again using Ginibre ensembles [1^, we obtained estimates of the first 500 mo- 
ments of the "balanced" variable |p^"^||p|. (The range of this variable is 



110592' 2562 J 

sec. 3.1].) In Figs. [3] and |4] we plot the ratios of these estimates to the exactly-known 
corresponding Hilbert-Schmidt moments [T3l App. D.6], given by 



(2n)!(l + «)2jl + 2a)2„ 



-n, a,a + |, —An — 1 



4-^3 



5a 



;1 



2i2n (3a + 1)^^ (6a + I) ^ \ -2n - a, -2n -2a,^-n 

(It is interesting to note, in this Hilbert-Schmidt setting, that the first moment of the 
balanced variable for the two-rebit systems [a = |], given by the formula is zero. For the 
two-qubit systems [a = 1], however, the analogous first moment is nonzero, that is, — 457^264') 
Further still, in Figs. |5] and |6] we plot the Bures-to-Hilbert-Schmidt ratios 
— — — ,n, = 0, 1, . . . , 24 of the bivariate moments. Let us note, interestingly, that 
the Hilbert-Schmidt bivariate moments themselves are symmetric/invariant under the inter- 
change of n and k, except when n or k equals zero [131 Apps. D.4, D.6]. The corresponding 
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FIG. 3: Ratio of numerically-generated Bures moments ((|/>^^||/9|)")^^^^*.^^-based upon 17,200,000 
random realizations-to corresponding exact Hilbert-Sclimidt moments. For n = 1, the HS moment 



IS zero. 
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FIG. 4: Ratio of numerically-generated Bures moments , ip,; /2-qubits 

random realizations-to corresponding exact Hilbert-Sclimidt moments 



)?^l!L«-based upon 17,200,000 



formula jTSj App. D.6] is 



where 



P \ \P\ 



HS 




HS 



(k + l)Jk + l + a)Jk + l + 2a)^ 

5a 



2n 



n, —k, a,a + |, —2k — 2n — 1 



—k — n — a, —k — n — 2a 



k+n fc+n — 1 ' 

2 ' 



;1 



HS 



k\ (a + l)j2a + l) 



(1) 



(3a+|)j6a + |),; 

In Fig.lTlwe display the (relatively well-behaved) ratio of the numerically-estimated Bures- 



ratio 



Bures 
2 — qubits 



Bures 
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to the HS-ratio 
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FIG. 5: Ratio of numerically-generated 6iuariate Bures moments ( p \p\ ) -based upon 

\ / 2—rebits 

10,500,000 random realizations-to corresponding exact Hilbert-Schmidt moments. For n = k = 1, 
the HS moment is zero. The HS moments themselves are symmetric under the interchange of n 
and k, except when n ox k equals zero. 

III. BURES DETERMINANTAL MOMENT COMPUTATATIONS 



A. Two-qubit results 

Now, let us present the following series \Pr)^^^^ubits' = 0, 1, ... 6) of ea;aci results 



(as inferred by us from high-precision computations). 



{\P'^\)T:1.U = -4 = " -0.00390625 (2) 
(the numerical estimate from the sample used to generate Fig. [2] being -0.003798959, while 



the Hilbert-Schmidt counterpart is {\p^'^\)^li = = — 22x3x^7x19 )' 



/I PT\ I \\Bures 137 137 , > 

\\P I \P\)2-quhits - ~ 1124597760 ~ ~2i9 x 3 x 5 x 11 x 13' ^' 



4576264 23x7x11x17x19x23^ 



/I PT\ I ,2\^''^^^^ 83 83 , , 

\\P I \P\ /2-gubits - ~ 4786288066560 ~ ~ 2^^ x3x5x7xllxl3xl9 ^ ' 

II pTl I ,^\Bures 565 565 , , 

\l^ I \P\ )2-qubits - "219011789856178176 ~ ~ 2^6 x 3 x 11 x 13 x 17 x 19 x 23 ^ ' 
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FIG. 6: Ratio of numerically-generated bivariate Bures moments ( \p \p\ ) -based upon 

\ / 2~qubits 

10,500,000 random realizations-to corresponding exact Hilbert-Schmidt moments. The HS mo- 
ments themselves are symmetric under the interchange of n and A;, except when n or k equals 
zero. 




FIG.7:R.foof<l-^^l!-&to'l^"l'-''>" 



\\P\ )2-rebit 



PT\ I :4\ Bures 



and 



139 



139 



2-qubUs 52562829565482762240 2^° x 3^ x 5 x 11 x 13 x 17 x 19 x 23 



(6) 



Bures 



14069 



14069 



2-qubits 1955073914750458962903040 252 x 5 x 13 x 17 x 19 x 23 x 29 x 31' 

(7) 
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^ X Bures 

Having also found (|p^"^| \p\ ) , A; = 6, . . . , 12, we were able to deduce (using the 

/ 2—qubits 



FindSequenceCommand of Mathematica) the underlying explanatory relation 

(8) 



, , Bures 

PT I I I 



1^ ' /2-quhits 8k'' + 36k^ - 82k^ - 68lP - 1366A; - 885 



^Bures ^38 (IG/c^ + 192k* + 883/^3 + 1947/c2 + 2062k + 840) 

2—qubits 



k{k{2k{2k{2k + 9) - 41) - 681) - 1366) - 885 
128(fc + 2)(fc + 2>){k + 4)(4fc + 5)(4A; + 7) ' 

fibres _ 315 {k + \)T{k+ I) r(2A; + 2) 



/2-g«bits V7rr(fc + 3)r(fc + 4)r (2A; + 1) ■ 



where 



The (lower-degree) Hilbert-Schmidt counterpart to the fundamental result ([8]) is [131 eq. 
(24)] 

PT I \k\ 

\P \\P\l^^-^_ k^ + 6P-k-A2 _ k{k{k + 6) - 1) - A2 



k\^^ 256A;3 + 3456fc2 + 15536A; + 23256 8(2/c + 9)(4A; + 17)(4A; + 19) ^ 



a=l 



(10) 

where (substituting a = 1 into ([T])) 

^k\^s 108972864000r(fc + l)T{k + 2)T{k + 3)T{k + 4) 



/..^ = r(4(fc + 4)) • 

In developing this last series of exact results (|2])-([8]), we combined symbolic methods with 
high-precision numerical computations as to the Bures expected values of monomials of the 
form 

At^A^^A^^'A^, (12) 

where the A's are the eigenvalues of p (with, of course, A4 = 1 — Ai — A2 — A3). If the 
monomial is fully symmetric (that is ki = k2 = k^ = k^) , then ^since |p| — AxA2A3A4^we can 
directly rely upon the determinantal moments implicit in the formulas [2|, [T^]. Otherwise, 
unfortunately-providing the greatest impediment to the research program pursued here-no 
general formulas seemed at hand. We initially, thus, had to proceed with the somewhat 
tedious process of numerically analyzing, to high precision, each specific monomial arising 
in our computations. For instance, in computing the two-qubit result ([S]) we have used (from 
0) {XiX2X3^A)2-qubits = lim = (2^ X 3 X ll)-i, as well as (our findings) 

(XlX2X3)T'Z. = —A>^l>^l)^''Z. = — , (13) 

\ i ^ •^/2—qubtts 50688 ^/2—qubits 16896 
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457 , ^.Bures 1127 

\^l^V2-aubits - 



'2-qubits 50688'^ ^/2-quhits 16896' 

However, we did succeed (again employing the Mathematica FindSequenceFunction com 
mand) in generalizing (sufficient for our needs in computing the explanatory formula ([8 
this last set of four results to 

{>^\^^>4^t^\~^)2-llbits {^i^2Mp\''~^)2-llhits ^k{k + 5) + 19 



l\„\k\Bures l\^\k\Bures Ak-2 ^ ^ 

h-auUts \\h'\ h-aubits ^'^ ^ 



(14) 



'2 — quhits \\r\ 1 2—qubits 

\k\Bures I \k \k \k \k\Bures 



where (|p|'=)r_",^:L. = {KA>^lK)rqubUs is given by Q . Further, 

{\\-^'\f'\t'\t')^!^i:ius _ k{k{Ak{k + 14) + 279) + 556) + 350 
{\pn2-7ulus k{k + 4){2k - l){2k + 1) 

{\'^^'\|\t'>^t')^^lSnts k{kmk + 19) + 439) + 996) + 770 



and 



(15) 
(16) 

(17) 



(|p|')f-;S.ts /c(/c + 4)(2fc-l)(2fc + l) 

{Xt^''A~^\t^\t^)^^luUts ^ {2k + 5){k{k{Ak{k + 32) + 1067) + 3148) + 2898) 
{\p\')2-i:Lts ~ k{k + A){2k-l){2k + l)^ 

B. Two-rebit results 

Turning from the generic (15-dimensional) two-bit case to the (9-dimensional) two-rebit 
one, we have found strong evidence that 

npr^.Bures ^ _ J663_ ^ _ 2663 ^ _o.00309594. (18) 

\/2~reb^ts 860160 2^3x3x5x7 ^ ^ 

The estimate for this quantity from the sample used to generate Fig. [T]is -0.00331101256, 
while {\p^'^\)^^i = —gig = — 2x3x11x13 ^PP- (It might be noted that a 



quite close approximation to the result (18)-but one we reject-is the reciprocal integer 
— 3I3 = —jY^jg ~ —0.00309597.) In computing this two-rebit result we have used (from [2]) 
{XiX2X3X4:)2-lebits ~ 8192 ~ 2"^^, as Well as our additional expected-value-monomial findings 



(cf. (13)) 



/ , 9 , , \ Bures 41 ,,o,o^ Bures 1399 ,,,, \ Bures 2507 ,,,, Bures 18463 

\ 1 2 ^ / 2-rebits 40960' ^ '^/2—rebits 286720' ^ / 2—rebits 286720' ^ ' 2—rebits 286720' 

(19) 



The computation to obtain ( 18 ) then took the form 



=-J^= (20) 

\/2-rebtts 860160 ^ ' 
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\ \Bures , / \ 2\ 2\Bures ?_/\3\ \Bures ^/\^\Bures , ^3 

-\^l^2^3 h-rebits onn \'^l'^2l 2-rebits ^ r \^1^2/2-refejts on 2-rebits 



\ 1 ^ J/ 2— reoits ' ggg \ 1 2/ 2— reoits \ 1 ^/ 2— reojts 2Q ^ 1 ' 2— rebjts ' 4Q9gQQ ' 

where the constant term = %{\i\2>^Ai)2Ztbits- 

Further, developing a two-rebit series analogous to the two-qubit series (|2|-([7]), we have 

/I PT\ I :\Bures 12T 127 /oi\ 

\|/^ I l/'l/2-refeiis - "1132462080 ~ ~ 223 X 33 X 5 ^ 

HS 



(while (|p^^||p|>::,=0) 



/ 1 PT I I 1 9 \ Bures 1257 3 X 419 

\\P I IPI /2-rebits ~ 1 nc^S97QQ/n T'i/l/in ~ 938 vC^v7vTT 



105827994173440 238 x 5 x 7 x 11 

py, 3,B„res _ 13813 __ 19x727 

(I/' I l/'l /2-rebits - 221661544159641600 ~ "2^7 x 3^ x 5^ x 7' ^ 



and 



I py, Bures _ 3927 _ 3x7x11x17 

I 1^1 /2~rebits - 28823037615171174400 ~ 2^0 x 5^ ' ^ 



Then, having found the next nine members of this series, we were able to generate the 
explanatory formula for the infinite series (cf. the two-qubit result (|8|, also the ratio of 
degree-5 polynomials in k) 

, \ Bures 
PT\ I ^\k\ 

(25) 



I /2-rebits _ 64A;5 + 128A;^ - 340A;3 - 1032A;2 - 1099A; - 384 

Bures 



where 



— " k (8/c2 -2k -1) (8A;2 + 18A; - 5) 

/ 2—rebits 

k{Ak{k{16k{k + 2) - 85) - 258) - 1099) - 384 
k{2k - l){2k + 5){4k - l){4k + 1) ' 

k\ Bures _ 322-8^r {2k + |) 

/2~rebits ~ v^(2P + 3A; + l)r(2A; + 4)' ^ ^ 



The (simpler, lower-degree) Hilbert-Schmidt counterpart to (25) is [131 eq. 3] 



, X HS 

P \P\ 



(A;-1)(A;(2A; + 11) + 16) 



32(A; + 3)(4A; + 11)(4A; + 13)' 



(27) 



It does appear that in the Bures case, such formulas are the ratios of degree-5n polynomials 
in k, while in the Hilbert-Schmidt case they are ratios of degree-3?2 polynomials in k |T3] . 
This " higher-degreeness" might be taken as some indication of the greater complexity of the 



Bures 2x2 determinantal moment problem (cf. (41), (42)). 
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As intermediate steps in deriving (25), we have developed the formulas 

/ \k+l \k \k \k—l\Bures /\2\ \ \ ^\k-l\Bures 

^2 ^3 ^4 /2-rebits _ \^1^2M\P\ l2-rebits 

l\„\k\Bures l\„\k\BuTes 

l2~rebits 1 2-rebits 

8A;2 + 22A; + 11 4 8 

+ 1, 



(28) 



(29) 
(30) 

(31) 



8F - 2A; - 1 4A; + 1 2A; - 1 
(At+^A^+^A^^A^^)f_";|g,,, _ Ak{k{lQk{k + 8) + 363) + 406) + 545 
(|P|')&. ~ {2k - l){2k + b){Ak - l){Ak + I) ' 

{X\'^^\lX't^\'l-^)^-lluts _ ^k{k{16k{k + 11) + 591) + 766) + 1325 
MWZ^s (2A;-l)(2A; + 5)(4fc-l)(4fc + l) ' 

and 

{Xi^^Xt'Xt^X'l-^)2-rtLts _ k{Ak{k{lQk{k + 20) + 1899) + 4678) + 20057) + 7680 
{\p\')^Zluts ~ /c(2/c-l)(2A; + 5)(4fc-l)(4A; + l) ' 

Further, 

, , 50654227 ^13 x 101 x 173 x 223 ^ 

^ /2-reUts 1307993702400 223 x 3^x52x7x11 ^ ' 

(The sample [Ginibre-ensemble] estimate-from Fig. [T]-was 0.000036949.) In attempting to 
find the n = 2 analogue of ( [25| , we have computed 

{X\XlXlXM')^Ztuts 8fc2 + 54fc + 87 



{W^^)^Zlms 8k' + 30k + 27- 



(33) 



(The right-hand side of this equation can be obtained by replacing k in (28) by k + 2.) Also, 

128A;5 + 2496A;^ + 19000A;3 + 70284P + 125922A; + 87213 



(Af AIA3 


P 


k\Bures 
1 2—rebits 


(|P 




2\Bures 
1 2—rebits 



128fc5 + 1344fc4 + 5080A;3 + 8700A;2 + 6642A; + 1701 
176 334 495 1260 5 ^ 



(34) 



63(4A; + 9) 6A; + 9 14A; + 7 44A; + 77 99(2A; + 9) 



We required fifteen distinct monomial-expected- value computations on the order of (33) and 



(34). Four of them could be directly obtained from (28)-(31) by the replacement of k by 
k + 2, and one immediately from the expected determinant formula, while the other ten 
required apparently new computations. We have been able to determine six of these ten so 
far. Among the successful computations were: 

{XiXlX3\p\'')^ZtMts _ 512A;6 + 13696^;^ + 143008A;^ + 751304^^ + 2096676^;^ -|- 2936118A; + 1596285 
{\p\''^^)2-r!bits ~ (4A; + 5) (128A;5 + 1344^^ + 5080A;3 + 8700^2 + 6642A; + 1701) 

(35) 
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5 675 _ 184 6480 14700 6 ^ 

858(2A; + 9) ^ 4A; + 2 ~ 3(4A; + 9) ^ 44A; + 77 " 52A; + 65 ~ 2A; + 3 ^ ' 

and 

/\5\2\ J^|fc\B 

-rebits 



/\5\2\ \^\k\Bures 
/2-rebits 



k+2\Bures 

(36) 

/2—rebi1 

512A;^ + 17792A;^ + 234144A;5 + 1576840^;^ + 5978260P + 12867318A;2 + 14589933A; + 6699810 



{k + 2){4k + 5) (128P + 1344^;^ + 5080A;3 + 8700^ + 6642A; + 1701) 
175 122500 17492 1936 55740 32765 1792 



+ 1. 



2574(2A; + 9) 39(4/e + 5) 63(4/e + 9) 6A; + 9 44/e + 77 84/c + 42 3(A; + 2) 
The Hilbert- Schmidt two-qubit counterpart of this last formula is 

+ 50k^ + 85 lA;^ + 6770A;=^ + 27234A;2 + 52970A; + 39084 
64 (1024A;6 + 30720^;^ + 383104^;^ + 2542080A;3 + 9465796A;2 + 18753960A; + 15444891) 

{k + 2){k + 3) {k{k{k{k + 45) + 620) + 3400) + 6514) 
~ 64(2A; + 9)(2A; + 11)(4A; + 17)(4A; + 19)(4A; + 21)(4A; + 23) ' 

The four formulas remaining to be determined are those for 

{\l\l\p\^)^Zttus (AfAi|p|'')f_— g,,, (AlA2|p|^)f— (A?|p|'=> 

l\n\k+'2\Bures 1 l\ „\k+'2\Bures 1 l\ r,\k+2\Bures aUxX 



Bures 



>2-rebits >2-rebits '2-rebit3 12-re.bits 

C. Computations for generic 3x3 real density matrices 

As a supplementary exercise in finding the expected values of monomials of eigenvalues 
for generic real n x n density matrices (p), we have been able to determine that for n — 3, 

(Ar^A^Ar^)f— ^ {xix2\pr')fz: ^ m^+sek+u 5 35 ^ 

|^|fc\^"^^* /iprt^"'"''^ 2 (8A;2 - 2/c - 1) 12A; + 3 6(2A; - 1) 

/ retrit \ / retrit 



and 

{Xi'^^\t'^3'')rZu _ 32P + 224fc2 + 308fc + 101 _ 4A:(8A;(A; + 7) + 77) + 101 

~ 32A;3 - 16A;2 - 2A; + 1 ~ (2A; - 1)(4A; - 1)(4A; + 1) 

/ retrit 



(38) 



(39) 



where 



/retrit " retrU ' r{3{k + l))r{k + 2)- ^^^^ 
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IV. DISCUSSION 



In our various Bures computations above, we have employed Euler-angle parameteri- 
zations of the two-rebit [HI App. I] and two-qubit density matrices pUj (while in the 
Hilbert-Schmidt study [13j, Cholesky-decomposition parameterizations were utilized). 

We note that in the three-dimensional simplex of eigenvalues (Aj > 0,S^^^Ai = 1), the 
Hilbert-Schmidt measure is proportional (a = ^ denoting the two-rebit systems, and a = 1, 
the two-qubit systems) to PI eq. (4.1)] 

PF^^U^.-^?''. (41) 



while [21 eq. (3.19)] 

n 



The normalization constants are, respectively, 80640 and 378378000 in the two HS cases of 
interest here, and ^ and for the two Bures cases, where the first member of each pair 
corresponds to a = ^ and the second member to a = 1. 

Further, it appears that the Hilbert-Schmidt determinantal moment formulas developed 
in [13] potentially provide certain "leverage" in obtaining (the desired) comparable Bures 
determinantal moment formulas. That is, if one had formulas for the expected values of 



the eigenvalue-monomials (12) in both the Bures and Hilbert-Schmidt cases, then one could 
reduce the computations for the Bures determinantal moments (avoiding or minimizing 
integrations over the [more numerous] non-eigenvalue parameters) by solving certain linear 
systems that incorporate the Hilbert-Schmidt counterparts. 

A referee has indicated that it would be of considerable interest to conduct similar anal- 
yses to those here for 3 x 3-dimensional systems, in particular in light of certain recent work 
pertaining to them. 
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